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ABSTRACT 

Aims. We present the results of a contemporaneous photometric and spectroscopic monitoring of two RS CVn binaries, namely 
A And and II Peg. The aim of this work is to investigate the behavior of surface inhomogeneities in the atmospheres of the 
active components of these systems which have nearly the same temperature but different gravity. 

Methods. The light curves and the modulation of the surface temperature, as recovered from line-depth ratios (LDRs) , are used 
to map the photospheric spots, while the Ha emission has been used as an indicator of chromospheric inhomogeneities. 
Results. The spot temperatures and sizes were derived from a spot model applied to the contemporaneous light and temperature 
curves. We find larger and cooler spots on II Peg (T sp ~ 3600 K) compared to A And (T sp ~ 3900 K); this could be the result of 
both the different gravity and the higher activity level of the former. Moreover, we find a clear anti-correlation between the Ha 
emission and the photospheric diagnostics (temperature and light curves). We have also detected a modulation of the intensity 
of the He I D3 line with the star rotation, suggesting the presence of surface features also in the upper chromosphere of these 
stars. A rough reconstruction of the 3D structure of their atmospheres has been also performed by applying a spot/plage model 
to the light and temperature curves and to the Ha flux modulation. In addition, a strong flare affecting the Ha, the He I D3, 
and the cores of Nal T>i t 2 lines has been observed on II Peg. 

Conclusions. The spot/plage configuration has been reconstructed in the visible component of A And and II Peg which have 
nearly the same temperature but very different gravity and rotation periods. A close spatial association of photospheric and 
chromospheric active regions, at the time of our observations, has been found in both stars. Larger and cooler spots have been 
found on II Peg, the system with the active component of higher gravity and higher activity level. The area ratio of plages to 
spots seems to decrease when the spots get bigger. Moreover, with the present and literature data, a correlation between the 
temperature difference AT = T p h — T sp and the surface gravity has been also suggested. 
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1. Introduction 

The rotational modulation of brightness and other pho- 
tospheric diagnostics in late-type stars is commonly at- 
tributed to starspots on their photospheres. The pres- 
ence of chromospheric inhomogeneities similar to the so- 
lar plages is pointed out by several stu dies based on the 
Ha an d /or other diagn o stics (s ee, e. g. . IStrassmeier et al. 
1993al ICatalano etall Il996l l2000l : iFrasca et all Il998l 
2000h . 



The simultaneous study of photospheric (spots) and 
chromospheric (plages) active regions (ARs) is important 
for a better understanding of the physical processes De- 
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curring during the emersion of magnetic flux tubes from 
below the photosphere. Indeed, the relative position and 
size of ARs at different atmospheric levels can provide 
information on the magnetic field topology. In the past 
ten years, rotational modulation due to surface inhomo- 
geneities at photospheric and chromosp heric level has been 
revealed in several RS CVn binaries. 



Catalano et al. 2000; Biazz o et al. 



|Frasca et al 
|2006Fk 



1998 



Simultaneous photometric and spectroscopic obser- 
vations have frequently shown a spatial a ssociation be- 
tween spots and plage s in RS CVn systems dRodono et al 



19871: ICatalano etall Il99rj IFrasca et all Il998l) as well 



as in young mildly-activ e solar-type stars ([Frasca et al 
2000$ Biazzo et al.l 2007 ). A close plage-spot association 



has been also detected in very active main-sequence 
single stars, like the r apidly rotating star LQ Hya 
(jStrassmeier et ah! 1993af). Also in close binary systems 
(e.g. TZ CrB. IFrasca et al.lll9971 ) and in extremely active 
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stars, like the comp onents of the contact binary VW Cep 
( Frasca et al . 1996) , there are evidences of spot-plage asso- 
ciations. Moreover, for some RS CVn binary systems there 
is some indication of a system atic longitude lag of 30°-50 ° 
between the plages an d spots (Catalano et al. ll996. 2000) 



curves have been published by Chugainovl ( 1976 ) Lon, 



truga 

term studies (e.g.. lHenry et al.lll995t iRodono et al.ll200 



is- 
50) 



In. recent works ([Catalano et al.l l2002t iFrasca et al 

20051 hereafter Paper I and Paper II, respectively), we 



presented the results of a contemporaneous spectroscopic 
and photometric monitoring of three active single-lined 
(SB1) RS CVn binaries (VY Ari, IM Peg, and HK Lac), 
showing that it is possible to recover fairly accurate spot 
temperature and size values by applying a spot model to 
contemporaneous line-depth ratios (LDRs) variation and 
light curves. Subsequently we have investigated the active 
region topology in t hese stars at both c hromospheric and 
photospheric layers (jBiazzo et al. 20061 hereafter Paper 
HI). 

In the present paper, a similar analysis is applied 
to two other active SB1 RS CVn's, namely A And and 
II Peg, which are stars with similar effective temperatures 
(T cff ~4700 and 4600 K, respectively ) but with different 
gravities, being \ogg ~ 2.5 for A And (Donati et al. 19951 ) 
and logg ~3.2 for II Peg (|Berdvugina et al.lll998bT ). 

A And (HD 222107) is a b right (V = 3j"82) and active 
giant, classified as G8 IV-III. ICalder (1938) was the first 
to discover its photometric variability, w hose amplitude 
sometimes reaches 0™30. Six years later, IWalkerl (| 19441 ) 
showed that A And is a SB1 with an almost circular or- 
bit of period 20^5212. It is an atypical member of the 
RS CVn class, because it is larg ely out of synchronism, it s 
rotational period being 531*952 (jStrassmeier et al.lll993bl ). 
Indeed, in most RS CVn binaries the rotational period 
of both components is equal to the orbital period of the 
system within a few percent. Thus A And, whose orbit 
is very close to be circul ar, is a puzzle for the theory of 
tidal friction (|Zahnlll977l ), which predicts that rotational 
synchronization in close binaries should precede orbit cir- 
cularization. It is one of the brightest of all known chro- 
mosph erically active binaries. Photoelectric ligh t curves of 
A And (|Bopp fe Noahlll98dlPoe fc Eatonlll985l ) show that 
its photometric period is somewhat variable, probably due 
to differential rotation and latitude drift of the spots dur- 
ing magnetic cycles. In addition, there is a long-term cycle 
of about 11 yr in the mean brightness. Phot ometric long- 
term studies have been performed also by iHenrv et al.1 
(1995). Moreover, the Ha emission has been found to be 
rotati onally modulated and anti -correlated with the light 
curve ( Baliunas fc Dupree 19821 ). 

II Peg (HP 2240 8 5) has been classified as a K2-3 IV- 
V SB1 by iRucinskil ( 1977 ) who noticed that its photo- 
me tric perio d was q uite close to 6?724183, as determined 
by lHallidavl ( 19521 ) for its essentially ci rcular orbit. It 
was classified as an RS CVn system by IVogtJ <ll981al fbh 
who found an amplitude of the V light curve of 0™43. 
II Peg is among the most active RS CVn binaries and be- 
longs to a small subset of binaries, including V711 Tau 
and UX Ari, in which H a always appears in emission 
( Nation fc Ramsevl 119811 ). The first photoelectric light 



have shown dramatic changes of the photometric wave, 
from almost sinusoidal, to irregular or flat. Moreover, 



iDonati et al.1 ( 1997t ) have clearly detected the magnetic 
field on this star. The so called flip-flop phenomenon (in 
which the dominant part of the spot activity changes the 
long itude every few years) has been re porte d in II Peg by, 
e.g., Berdvugina fc Tuominen I ( 1998) and Rodono et al 



(2000) and theoretically analyzed bv lElstner fc Korhonen 
(120051) . 



The aim of the present work is to investigate the 
starspot characteristics of this two active stars with the 
same temperature but different gravity and activity level, 
as well as to study the location of the excess Ha emission 
and the degree of spatial superposition between surface 
inhomogeneities at different atmospheric levels. 

2. Observations and reduction 

2.1. Spectroscopy 

Spectroscopic observations have been performed in 1999 
and 2000 at the M. G. Fracastoro station (Serra La 
Nave, Mt. Etna) of Catania Astrophysical Observatory 
with FRESCO (Fiber-optic Reosc Echelle Spectrograph 
of Catania Observatory), the echelle spectrograph con- 
nected to the 91-cm telescope through an optical fiber 
with a 200-/zm core diameter. The spectral resolution was 
R = A/AA ~ 14000, with a 2.6-pixel sampling. 

The data reduction was performed by using the 
echelle task of the IRAI0 package following the stan- 
dard steps of background subtraction, division by a flat 
field spectrum given by a halogen lamp, wavelength cali- 
bration using the emission lines of a Thorium-Argon lamp, 
and normalization to the continuum through a polynomial 
fit. 

We have removed the telluric water vapor lines at 
the Ha wavelengths using the spectra of Altair (A7 V, 
vsvcii = 245 kms -1 ) acquired during the observing runs. 
These spectra have been normalized, also inside the very 
broad Ha profile, to provide valuable templates for the wa- 
ter vapor lines. An interactive procedure, allowing the in- 
tensity of the template lines to vary (leaving the line ratios 
unchanged) until a satisfactory agreement with each ob- 
served spectrum is reached, has been applied to correct the 



obser ved spectra for telluric absorption (see IFrasca et al 
120001) . 

2.2. Photometry 

The photometric observations were performed in the B, V, 
and R Johnson filters at the Ege University Observatory. 
The observations were made with an un-refrigerated 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under 
cooperative agreement with the National Science Foundation. 



A. Frasca et al.: Spots, plages, and flares in A And and II Peg. 



3 



Hamamatsu R4457 photometer attached to the 48-cm 
Cassegrain telescope. 

A And was observed from July 9 to November 2, 1999 
for a total of 35 nights, using ip And and k And as com- 
parison (C) and check (Ck) star, respectively. II Peg was 
observed from July 3 to November 21, 2000 for a total of 
25 nights, using HD 224083 and BD+27°4648 as compar- 
ison and check star, respectively. 

The differential magnitudes, in the sense of variable 
minus comparison (V— C), were corrected for atmospheric 
extinction using the seasonal average coefficients for the 
Ege University Observatory. The light curves were ob- 
tained by averaging individual data points taken in the 
same night (from 5 to 20). The standard deviation of each 
observed point, as measured from the differential magni- 
tude V - C and C - Ck, ranges from ±(F005 to ±0 m 015. 

3. Data analysis 

3.1. Ha line analysis 

The hydrogen Ha line is one of the most useful and easily 
accessible indicators of chromospheric activity in the op- 
tical spectrum. Furthermore, it is very effective, both in 
the Sun and in active stars, for detecting chromospheric 
plages, due to their high contrast against the surrounding 
chromosphere. 

II Peg is one of the few RS CVn stars which dis- 
plays the Ha line always in emission above the local 
continuum. For the majority of the active binaries, in- 
stead, only a filling-in of the Ha core can be seen. In 
these cases, the "spectral synthesis" method, based on the 
comparison with synthetic spectra from radiative equilib- 
rium models or observed spectra of non-active standard 
star s (reference spectra), has been successfully used (see 



a possible contamination by chromospheric emission in 
the Ha core of the reference stars, could be still present. 
However, this contribution can be completely neglected 
for very active stars, like the RS CVn systems here in- 
vestigated. More inform ation about the spectral sy nthesis 
method can be found in 



,tion apout the s pectral sy: 
Frasca fe Catalanol (|l994[ ) 



3.2. The Helium D3 line 

As a further indicator of chromospheric activity we have 
analyzed also the Hei A 5876 (D3) line, which is nor- 
mally seen as an absorption feature in active stars (e.g., 
lHuenemoerderlll98d iBiazzo et al.ll2006ll2007t ). The refer- 
ence spectra adopted do not show any He 1 absorption nor 
emission, as expected for non-active stars. In this case, 
the spectral subtraction technique allows to emphasize 
the He 1 line, cleaning the spectrum from nearby photo- 
spheric absorption lines, and to measure its equivalent 
width, EWsei ■ For EWrcI we have adopted the usual con- 
vention that an absorption line has a positive EW . 

The helium A 5876 line, because of its high excitation 
potential, is a good tracer of the regions of higher temper- 
ature and excitation in solar and stellar chromospheres. 
Recent models seem to indicate that the primary mech- 
anism responsible for the formation of the He 1 triplet is 
the collisional excitation and ionizati on by electron impact 
follow ed by recombination cascade (jLanzafame fc Byrne 
1995ft . 

In the Sun, the He 1 D3 line appears as an absorption 
feature in p lages and w eak flares and in emission in strong 
flares (e.g., IZirinl ll988L In active stars, the Hei D3 is usu- 
ally observed in absorption and sometimes in emission, 
like during fla r e events in the most active RS CV n star s 



(jMontes et all Il997l Il999t iGarcia- Alvarez et all l2003h . 



( reference spectra), h as be en successfully used (see , " " " 1 1 1 """ " ^ "*i 

H^rbkl[l985llBardenl EUl iFrasca & Cataknol fl99l This is related to the electronic temperature and density in 



MontesetaL 1995). The difference between observed and 



reference spectrum provides, as residual, the net chromo- 
spheric Ha emission, which can be integrated to find the 
total radiative losses in the line. 

Since II Peg and A And are SB1 systems, we used only 
one standard star spectrum to reproduce their observed 
spectra. This standard star spectrum has been rotation- 
ally broadened by convolution with a rotational profile 
with vsini = 22.6 kms -1 for II Peg and wsini = 6.5 



the e mitting region (e.g., IZirin 1988; lLanzafame fc Bvrnd 
19951 ) . A contribution to the emissivity from overionisa- 



kmf 



for A And to mimic the active star in absence of 



chromospheric activity. Figures Q] and [3] show samples of 
spectra in the Ha region of A And and II Peg, respec- 
tively. The residual Ha equivalent width, AEWna, has 
been measured by integrating all the emission profile in 
the difference spectrum (see the lower panel of Fig. [T]). 

The error, a(AEWn a ), has been evaluated by mul- 
tiplying the integration range by the photometric error 
on each point. This latter has been estimated by the stan- 
dard deviation of the observed flux values on the difference 
spectra in two spectral regions near the Ha line. Although 
the errors a(AEW-n a ) are generally a few hundredths of 
an A, a very small systematic error in AEWb_ u , due to 



tion due to coronal EUV back-radiation is expected to 
play a role when the transition region pressure is below 1 
dyne cm -2 ( Lanzafame fc Bvrneill995h . This contribution 
is therefore relevant for moderately active stars. In any 
case, our analysis in independent on the detailed mecha- 
nism of formation. 



3.3. Photospheric temperature from LDRs 

Line-depth ratios (LDRs) can be used for detecting the 
temperature rotational modulation in active RS CVn 
stars. Such diagnostics allow to detect temperature vari- 
ations as small as 10-20 K at the resolution of our spec- 
tra and with a good signal-to-noise ratio (S/N > 100). 
The precision of this method is improved by averaging 
the results from several line pairs, as discussed in Papers I 
and II. Typically we used from seven to fifteen line pairs, 
depending on the star's vsini, in the wavelength region 
around 6250 A to produce an average value of the star 
photospheric temperature at each phase. 
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Fig. 1. Top panel: Example of observed spectrum of A And 
at the rotational phase (f> = 0?761 in the Ha region (full 
line) together with the non-active template (dotted line). 
Lower panel: Difference between observed and template 
spectra emphasizing the core chromospheric emission. 

We found a clear rotational modulation of the photo- 
spheric temperature recovered from LDRs for three SB1 
RS CVn binaries, namely VY Ari, IM Peg, and HK Lac 
(Paper I). In Paper II we showed that these temperature 
curves are very well correlated with the contemporaneous 
light curves. Moreover, the simultaneous modeling of the 
temperature and light curves enabled us to derive temper- 
ature and size of the starspots. 

In principle, rotational line-broadening must be taken 
into account in a technique such as the LDR. However, 
v sin i is usually smaller than the instrumental resolution 
and rotational broadening can be safely ignored. This 
is not the case for II Peg, whose rotational broadening 
(wsini = 22.6 kms -1 ) is larger than the resolution of 
our spectra. In this case we have therefore produced a 
LDR-T e ff calibration using rotationally broadened refer- 
ence spectra. The values of temperature deduced for II Peg 
with these calibrations are, however, close to those found 
with the calibration at vsini = reported in Paper I. 

Further information about the LDR method can be 
found in Papers I and II. 

4. Results 

4.1. A And 

For A And, nine LDRs have been used and transformed 
in temperature variations by means of the LDR-T e g cali- 
bration with no rotational broadenin g. The rotational ve - 



locity of A And, namely 6.5 kms (jDonati et al.l 119951 ), 



kms^ 1 , so that no correction for rotational broadening 
is needed. All the LDRs converted into temperature and 
combined in a single temperature curve lead to a fairly 
well-defined temperature variation as a function of the ro- 
tational phase and well correlated with the optical light 
curve (Fig. [2]). Unfortunately, the very long rotation pe- 
riod has prevented us to obtain a complete phase coverage, 
but we have enough data around the maximum and min- 
imum of the curve for performing a meaningful analysis. 

The rotational phases have been derived according to 
the following ephemeris 

HJD^o = 2 443 829.2 + 53 d 952 x E , (1) 

where the i nitial heliocentric Julian day is that of 
Henry et al. (Il995) an d the r otational period is taken 



from Strassmeier et al. (Il993ri . The temperature maxi- 
mum, with a value of 4740 K, occurs at phase </> ~ 0?65. 
The full amplitude variation of the effective temperature 
is AT e ff ~ 80 K, corresponding to about 2% of the average 
temperature. The light curve displays a variation ampli- 
tude AV = m 225. 

fj, Peg (G8 III), with B - V = 0'. n 934, has been used as 
non-active Ha reference, because it is of the same spectral 
type and nearly the same color as A And. The Ha line of 
A And is always in absorption but with a variable filling-in 
in the core (see Fig. Q] for an example) . The values of the 
residual emission AEWn a with their errors are plotted in 
Fig. [2] as a function of the rotational phase. Two AEWn a 
values taken in the same nights (marked with an arrow) 
could be related to a mild flare event. In these hypothesis, 
we can only give an upper limit of two days for the flare 
duration, because there is no sign of Ha enhancements ei- 
ther in the night before or in that following the event. It is 
also very difficult to estimate the flare energetics because 
we don't know the flare duration and the Ha level at the 
flare peak. However, we have roughly evaluated the energy 
released in the Ha line during the event assuming a peak 
value AEWn a — 0.37 A (underestimate) and a flare dura- 
tion of 2 days (overestimate). We have converted the net 
Ha equivalent into luminosity by means of the equation: 

Ln a = L^AEWn a = 47rd 2 f 6563 AEW Ra 



= And 



Ft 



(i5(i:S 



1 (-°.4Vo-8.451) AWhq) 



(2) 



-F5556 

where ^6563 an d /6563 arc the luminosity and the flux at 
Earth of the continuum at A = 6563 A, d is the distance 
and io(-°- 4y o-8.45i) i s the Earth flux at A = 6563 A from a 
star of Vo (de-reddened) magnitude. The continuum flux- 
ratio 



-F5556 I ; " II -1 

low-resolution spectra (jHauschildt et al.H1999l ). We find a 
luminosity at peak Lu a ~ 3 x 10 30 ergs -1 and an approx- 
imate value for the total energy emitted in the Ha line of 
E Ka « 8 x 10 34 erg. 

Besides this possible flare event, a clear anti- 
correlation between temperature and Ha emission is ap- 
parent with a similar shape of the curves. Hei D 3 ab- 
sorption has been also detected, but no clear modulation 
emerges from the scatter of the data. 



has been e valuated using Next Ge n synthetic 



is in fact lower than the FRESCO resolution of about 7 
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A And 




Table 1. Temperature values and parameters of the sub- 
tracted spectra of A And. 





= 0.10 



= 0.35 



Fig. 2. From top to bottom. Johnson V photometry, T e g, 
AEWiia , and EWhgI as a function of the rotational phase 
for A And. The synthetic temperature, light, and EWn a 
curves, calculated with our spot/plage model (see Sect.E]), 
are reproduced with full lines in each box. The inset dis- 
plays the AEWn a versus Julian Day. A schematic map of 
the spot and plage distributions, as seen at two different 
rotational phases, is also shown in the bottom. 



The values of T off , AEWn a , and EWnei for all the 
observed spectra are listed in Table [1] 

4.2. II Pegasi 

Seven line-depth ratios have been used for II Peg. A few 
line pairs have been disregarded due to severe blending 
at the tisini of II Peg. The LDR variations have been 
transformed into temperature variations by means of the 
calibration based on standard stars rotationally broad- 
ened at the same y sin % of II Peg, namely 22.6 kms -1 
(Ber dvugina et al.lll998a ). The rotational modulation of 



the temperature and V magnitude is shown in the two 
upper boxes of Fig. [4] as a function of the phase calculated 
according to the following ephemeris 



H JD 


Phase 








(+2 400 000) 




( K ) 


(A) 


(A) 


51 379.527 


0.945 


4683± 7 


0.22±0.03 


0.063±0.018 


51 379.535 


0.945 


4676± 7 


0.21±0.03 


0.074±0.030 


51 380.500 


0.963 


4675± 6 


0.25±0.03 


0.068±0.034 


51 380.512 


0.964 


4672±21 


0.24±0.05 


0.089±0.031 


51 381.527 


0.982 


4665± 8 


0.29±0.03 


0.092±0.043 


51 381.539 


0.983 


4678±13 


0.28±0.02 


0.059±0.018 


51 382.531 


0.001 


4661± 4 


0.28±0.02 


0.056±0.017 


51 382.539 


0.001 


4660± 8 


0.27±0.02 


0.049±0.015 


51 384.480 


0.037 


4665±14 


0.28±0.04 


0.053±0.012 


51 385.527 


0.056 


4689± 7 


0.30±0.03 


0.066±0.029 


51385.535 


0.057 


4687± 6 


0.30±0.02 


0.083±0.027 


51386.512 


0.075 


4701± 9 


0.27±0.02 


0.066±0.024 


51386.520 


0.075 


4697±10 


0.27±0.03 


0.080±0.025 


51387.520 


0.093 


4684± 6 


0.38±0.03 


0.073±0.024 


51387.527 


0.094 


4683± 5 


0.36±0.03 


0.057±0.018 


51388.523 


0.112 


4672± 3 


0.27±0.04 


0.078±0.021 


51388.531 


0.112 


4686± 5 


0.26±0.03 


0.083±0.021 


51389.488 


0.130 


4686± 7 


0.28±0.03 


0.087±0.021 


51389.496 


0.130 


4677± 6 


0.29±0.04 


0.089±0.022 


51390.488 


0.148 


4685± 6 


0.25±0.04 


0.072±0.020 


51390.500 


0.149 


4700± 9 


0.24±0.06 


0.077±0.019 


51407.586 


0.465 


4736± 7 


0.11±0.02 


0.041±0.018 


51417.520 


0.650 


4729±11 


0.15±0.03 


0.064±0.023 


51417.590 


0.651 


4738± 4 


0.13±0.03 


0.062±0.023 


51421.578 


0.725 


4719± 5 


0.18±0.05 


0.070±0.022 


51423.520 


0.761 


4735± 5 


0.13±0.02 


0.063±0.020 


51426.547 


0.817 


4716±12 


0.16±0.03 


0.075±0.021 


51426.555 


0.817 


4711± 9 


0.17±0.04 


0.067±0.020 



wher e the rotational pe r iod and the initial epoch are take n 
from lBerdvugina et alj (|l998al ) and lHall fe Henrvl (| 19831 ). 
respectively. The temperature variation is about 3%, with 
an amplitude AT e ff ~ 130 K, while the light curve am- 
plitude is AV ~ 0^63. The two curves seem to have a 
slightly different shape and a small shift of the phase of 
minimum. 

S Eri (KO IV) has been used as a reference inactive star, 
whose spectrum has been broadened at the rotational ve- 
locity of II Peg. A sample of II Peg spectra in the Ha re- 
gion at different phases is shown in Fig. [3l The Ha line is 
always in emission above the continuum level with a small 
core reversal from <j> = 0?157 to 0?304, while a pure emis- 
sion with asymmetric shape is observed at other phases. A 
very strong and broad emission is observed at <f> = 0?906. 
The AEWhq and EWaei as a function of the rotational 
phase are plotted in Fig. [4] and listed in Table [2] together 
with the temperature values. In Fig. [4j a sudden increase 
of the Ha equivalent width at HJD=2 451 865 (0 = 0?906), 
which corresponds to the onset of a flare as witnessed 
also by the Ha profile, is marked with an arrow. This 
behavior is also present in the equivalent width of the He I 
line, which at that epoch was in emission (Fig. [5|), unlike 
outside-flare spectra in which it is always in absorption. 
The strength of the event is also witnessed by the remark- 
able filling-in of the sodium Di and D2 lines (Fig.[5l), sim- 
i larly to what found for HR 1099 by[G arcia- Alvarez et al 



HJD 



4>=o 



2 443 030.24 + 6? 724333 x E, 



(3) 



(|2003l ) . Since II Peg was in a quiescent phase on November 
15 th , 2000, it is possible to estimate a lower limit of two 
days for the life of the flare, assuming that the outburst 
was composed by only one event. In fact, one day after 
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6555 



6570 



Fig. 3. Time series of the Ha profile of II Peg (full lines) 
together with the non-active template (dotted line). Note 
the very strong Ha profile with broad wings observed dur- 
ing a strong flare on November 16 th followed by a nearly 
"normal" emission profile the night after (thin line). 



the flare peak (November 17 th , HJD=2 451866) there is 
still some Ha emission excess, with respect to the average 
modulation curve (big circle in Fig. |4|) . Unfortunately, our 
photometric data does not include the [/-band, sensitive 
to flares, and therefore no comparison with photometry 
can be made. 

We evaluated the energy released in Ha according to 
Eq. [21 Integrating the "excess luminosity" on all the du- 
ration of the flare ( ^ 2 days) , the total energy emitted in 
the Ha line is ^ 10 35 erg. It is worth noticing that 
this flare has occurred near the V minimum, which de- 
notes a possible spatial association with the photospheric 
spots, in analogy with the most energetic solar flares, 
the so-called "two-ribbon" flares occurring in the biggest 
and more complex spotted areas. Analogously to the Ha 
line, we estimated the peak luminosity for the He I line 
as £h~c ~ 9 x 10 28 ergs -1 and the total emitted energy 
as E]£* < 10 34 erg. It seems that the flare had a shorter 
duration in the He I line, because the £Wh~cI of the last 
spectrum is at the same level as in the quiescent phase. 






= 0.00 



= 0.70 



Fig. 4. From top to bottom. Johnson V magnitude, aver- 
aged T e ff, total Ha emission and Hei D3 equivalent width 
as a function of the rotational phase for II Peg. The syn- 
thetic temperature, light, and EWho curves, calculated 
with our spot /plage model (see Sect. [5]), are reproduced 
with full lines in each box. The inset displays the AEWn a 
versus Julian Day with an arrow marking the flare peak. 
A schematic map of the spot and plage distributions, as 
seen at two different rotational phases, is also shown in 
the bottom. 



5. The spot/plage model 

We have shown in Paper II that, with a spot model ap- 
plied to contemporaneous light and temperature curves, 
it is possible to reconstruct the starspots distribution and 
to remove the degeneracy of solutions regarding spot tem- 
perature and areas. Our spot model assumes two circu- 
lar active regions whose flux contrast (-F S p/-Fph) can be 
evaluated through the Pl anck spectral energy distribu- 
tion (SEP ) . the ATLAS9 (|Kurucd ll993l) and PHOENIX 



Next Gen (Hauschildt et al 



1999) atmosphere models. 



In Paper II we have evaluated temperature and sizes 
for the starspots observed on VY Ari, IM Peg, and HK Lac 
in the fall of 2000. We have also shown that both the at- 
mospheric models (ATLAS9 and NextGen) provide values 
of the spot temperature, T sp , and area coverage in close 
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Table 2. Temperature values and parameters of the sub- 
tracted spectra of II Peg. 



HJD 


Phase 


T off 


AEW Ha 


BWhci 


(+2 400 000) 




(K) 


(A) 


(A) 


51856.359 


0.564 


4546± 4 


2.54±0.11 


0.012±0.015 


51858.328 


0.857 


4478±29 


3.15±0.21 


0.028±0.038 


51 859.422 


0.020 


4525±34 


2.83±0.13 


0.019±0.021 


51860.348 


0.157 


4525±11 


2.24±0.12 


0.032±0.028 


51861.328 


0.303 


4573± 7 


1.84±0.11 


0.030±0.017 


51862.379 


0.460 


4570±25 


2.10±0.12 


0.034±0.018 


51863.383 


0.610 


4514±20 


2.62±0.12 


0.047±0.025 


51 864.488 


0.773 


4446±20 


3.05±0.13 


0.013±0.038 


51865.379 


0.906 


4453±10 


5.19±0.15 


-0.167±0.051 


51866.457 


0.066 


4495±23 


3.00±0.16 


0.020±0.027 



1.2- 




5870 5880 5890 5900 5910 

A (A) 



Fig. 5. Two observed spectra of II Peg in the region of 
He I D3 and Nai Di,D2 at the flare peak (full line) and 
in a quiescent phase (dotted line). Note the strong Hei 
emission and the remarkable filling-in of the Nai doublet 
during the flare. 



Table 3. Spot/plage configuration for A And a . 



Radius 


Lon. 


Lat. 


Tap/T ph 




A rel 


(°) 


n 


(°) 


(K) 




Spots 


26.5 


343 


57 


815+°- UB4 


3861±™* 


°- 076 -0.014 


17.1 


64 


9 








Plages 


22.5 


360 


63 






0.068 


19.8 


48 


18 









a Limb-darkening coefficients /^v— 0.797 and /X6200— 0.68, T p h — 4738 K, 
A_EW quiot =0.114 A, and 

-^plagc 

/-Fchrom = 8 have been used. 



such a low flux ratio and they could not reproduce the 
observed modulations. On the other hand, very high val- 
ues of flux ratio (-Fpiage/^chrom ~ 10) would imply very 
small plages producing top- flattened modulations that are 
not observed. So the flux ratio was fixed to values in the 
range 3-8, that is also typical of the bri ghtest parts of so- 
lar pla ges or of flare regions. Note that lLanzafame et al. 
(j2000h find F plago /F chrom ~ 3 for HR 1099, which has T eff 
and logg similar to II Peg. 

The solutions essentially provide the longitude of the 
plages, giving only rough estimates of their latitude and 
size. We have searched for the best solution by varying the 
longitudes, latitudes and radii of the active regions. The 
radii are, however, strongly dependent on the assumed 
flux contrast i^piage/^chrom- Thus, only the combined in- 
formation between plage dimensions and flux contrast, i.e. 
some kind of plage "luminosity" in units of the quiet chro- 
mosphere (ipi a g e /L qu i c t) can be deduced as a meaningful 
parameter. Note also that we cannot estimate the true 
quiet chromospheric contribution (network) , since the Ha 
minimum value, AEW qu i c t, could be still affected by a 
homogeneous distribution of smaller plages. 



agreement, while the black-body assumption for the SED 
leads to underestimate the spot temperatur e. This result 



is also in agreement with the findings of lAmado et al 
(1999). Since we have no long-term record of the pho- 



tosphcric temperature, we have assumed the maximum 
value obtained in our observing runs as the "unspotted" 
temperature for the modeling. 

In order to analyze the chromospheric rotational mod- 
ulations, we have extended our spot model, allowing also 
for bright active regions, with the aim of further inves- 
tigating the degree of spot-plage correlation. Given the 
scatter in the data, two bright spots (plages) are fully suf- 
ficient to reproduce the observed variations. 

A parameter that must be fixed in the model is the 
flux ratio between plages and surrounding chromosphere 
(-Fpiage/^chrom)- Values of F p i agc /F chrom « 2, that can be 
deduced averagin g solar plages in Ha (e.g., ISvestkalll976t 
lAyres et al.l ll98(j). are too low for modeling the high am- 
plitudes of Ha emission curves observed for A And and 
II Peg. In fact, extremely large plages, covering a signifi- 
cant fraction of the stellar surface, would be required with 



5.1. X Andromedae 

From the unspotted temperature and magnitude of A And, 
namely T p h = 4740 K and the historical light maximum 



Vmax = 3'. n 70 (jBovd et al .1119831 ). and the Hipparcos par- 



allax of 38.74 ±0.68 mas, the derived star radius was 
R = 7.51 Rq. The inclination of the rotation axis with 
respect to the line of sight de rived from this valu e of star 
radius, the v sin i=6.5 km s _1 (jDonati et al ]| 19951). and the 



rotati on period comes out to be i ~ 67°. [Donati et al 
([19951 ) found the same value of the star radius and esti- 
mated a n inclination i = 60 °li°. Since our i value is the 
same as ( Donati et al.lll995T) within the errors, we adopted 
i = 60° for the spot modeling. 

The results of the grids of solutions for the light 
curve and the temperature curve are displayed in Fig. O 
while the spot/plage configuration is given in Table [3] 
and displayed in Fig. [2] which also shows the synthetic 
curves superimposed to the data as full lines. In this case 
Fpiage/Fchrom was set to 8. A smaller contrast value would 
give rise to very big plages with a worse fitting of the Ha 
curve. 
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A And 
♦ 



0.20 : 
0.15 : 
0.10 : 
0.05 : ' 



P ♦ 



• « • 



0.5 0.6 



0.8 0.9 



1 .0 



Tsp/Tph 



Fig. 6. Grids of solutions for A And. The filled circles rep- 
resent the values of spot temperature and area for the 
light-curve solutions, while the diamonds represent the so- 
lutions for the temperature curve. The hatched area, in 
each box, is the locus of the allowed solutions accounting 
for the data uncertainty. 



For A And, IBopp fc Noahl (|l98dh and IPoe fe Eaton 
(1985) found that the asymmetric shape of the light curve 
requires two spots at different longitudes and that these 
spots are 800-1050 K cooler than the surrounding pho- 
tosphere, in close agreement wi th our findings. The two 
spots, revealed in the images of iDonati et al.l (1995) and 
covering altogether about 12% of the total stellar surface, 
had temperatures of 4000±300K, while the photospheric 
temperature was determined to be 4800±50K, giving a 
temperature difference AT ~ 800 K. They also calculated 
that the strong global magnetic fields recently detected 
in RS CVn systems could provide strong enough mag- 
netic braking to explain the o bserved non-synchron i zation 
of A And's rotation. Finally. [Padmakar k Pandev (1999) 
find AT = 800 ± 30 K and A rel = 8.5 ± 0.3 % setting 
T ph = 4800 K and i = 60°. The value of AT ~ 850 K 
derived by us is in close agreement with all these previ- 
ous determinations, but lO'Neal et al. (1998b) find cooler 
spots (T sp ~ 3650 K) by using the TiO bands. 

5.2. II Pegasi 

For II Peg we adopted unspo tted magnitude and tem- 
perature values of Vm ax = 6™9 iBerdvugina et al. I (|l998bh 
and r p h = 4573 K (present work). From the unspotted 
and de-reddened V magnitude and the Hipparcos paral- 
lax of 23.62 mas, a star radius of 2.76 Rq is derived with 
the Barnes- Evans relation. From the radius, the rotational 
period of 6"?724, and the vs'mi value of 22.6 kms -1 , an 
inclination of 60°!iq has been deduced. The mass we de- 



rived from the evolutionary tracks of Girardi et al.l (12 0001 



for the primary component is 0.9 M®. IBerdvugina et al 



(|1998af) derived a mass of 0.8 ± 0.1 M Q and a radius of 
3.4±0.2 Rq (log g ~ 3.2) for the primary component and a 
mass of 0.4±0.1 Mq for the unseen secondary component, 
whose spectral type has been estimated to be M0-3V. 



The result of the intersection of the grids of solutions 
for the light curve and the temperature curve (Fig. [7]) 
provides a T sp /T p h of 0.787, i.e a T sp ~ 3600 K. However, 
the spot model with NextGen SED has been also applied 
and similar values to those presented in Table [4] have been 
found. 

To reproduce the AEWn a curve a relative flux 
-Fpiagc/^chrom > 2 was needed. A value of F p i age /F c hroin = 
4 gives a satisfactory fit of the data, though, as previ- 
ously outlined, this parameter cannot be constrained by 
the modeli ng of the Ha curv e. H owever, the analysis car - 
ried out bv lBusa et al. ( 19991 ) and Lanzafame et al. ( 2000h 
on HR 1099, which has T e g and log g similar to II Peg, 
points to a contrast of Fpiage/Fchrom ~ 3. Such analysis, 
based on the modeling of the Mgn h k k and Ha profiles, 
gives more constraints on the the flux contrast and the 
AR area. 

The configuration of the active regions is represented 
in Fig. 0] and the active region parameters are listed in 
Table gj 

The first work on s pot modeling for II Peg was 
made by IBopp fe Noahl (|l980h . who showed that, for a 
satisfactory modeling of their asy mmetr ic light curve, 



two cool spots were s u fficien t. IVogtl (|1981al) and 



Huenemoerder fc Ramsevl (|1987l ) made a quantitative 



study of the effect of spots in the TiO bands and 
found that a substantial fraction of the photosphere 
must be spo tted (with a coverage of 35-40%). Then, 
VogH d!981bh found that the light and color curves 



could be reproduced with a cool spot having an ef- 
fective temperature of 3400±100 K and covering about 
37% of one hemisphere on the star. In subsequent years, 
several other authors have produced spot models for 
several datasets obtained in different epochs and de- 
rived spot temper atures. Among the s e, the more relevant 
works a re those of Nation fc Ramsey (Il981 ^ Poe fc Eaton 
(Il985l). iRodono et all (|l986l ). iBvrne fc Marand (jl987l ). 
an dlBovd et al.l (Il987l) . 



Donati et al.l (|1997f ). by means of the Zeeman-Doppler 



Imaging technique, which takes advantage of the Dopplcr 
effect in rapidly rotating stars to separate in wavelength 
the disk-integrated polarized profiles of magnetic lines, 
detected clear signatures of magnetic field with possi- 
ble conce ntrations near the longitudes of the spots ob- 
served by Henry et al. ( 19951 ). Moreover, it was shown for 
this star that the changes in the light curve of II Peg 
are consistent with re-arrange ments of th e spot distri- 
bution over the stellar surface. Neff et al. ( 19951 ). from 
TiO molecular bands, found that cool starspots (T sp = 
3500±200 K) are always visible, with a fractional pro- 
jected coverage of the visible hemisphere varying f rom 
54% to 64% as the star rotates. lO'Neal fc Nefll <jl997l) de- 
tected excess of OH absorption due to cool spots on the 
surface of II Peg and found a spot filling factor of 35- 
48 %, consistent with t h e minimum va lue of 40% found 
bv iMarino et all (|l999h . iHatzed (|l995l ) revealed polar or 
high-latitude spots and severa l equatorial s p ots wi th a 
total coverage of about 15%. lO'Neal et al.l (|1998al ) re- 
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ported the first spectroscopic evidence for the multiple 
spot temperature for this star. From TiO molecular band 
observations, they found that spot temperature varied be- 
tween 3350 to 3550 K during the epoch of September 1996 
to October 1996, whereas the starspot filling factor was 
constant (about 55%). Finally, Berdvugina et alJ 1 1998bl ) 
found, from their surface images, that the high-latitude 
spots were the m ajor contribution to the ph otospheric 
activity of II Peg. Padmakar fc Pandev ( 19991 ) obtained 
AT = 740 ± 45 K an d = 15.7 set ting T ph = 4300 K 
and % = 34°. Recently. IGu et alJ (|2003l ). from the Doppler 
imaging analysis, have found for II Peg changes in spot dis- 
tribution, including the position, intensity and size of the 
spots occurring in 1999-2001. In particular, they find in 
February 2000 a large high-latitude spot and a small weak 
low-latitude spot, i.e. the same spot configuration found 
in the present work for the observing season of November 
2000. 

Table 4. Spot/plage configuration for II Peg a . 



Radius 


Lon. 


Lat. 








(°) 


n 


(°) 


(K) 




Spots 


47.8 


320 


60 


U-787_ 041 


3599+ K 


0.189_ 018 


16.8 


202 


8 








Plages 


35.9 


322 


60 






0.132 


18.8 


206 


8 









a Limb-darkening coefficients /j.v— 0.836 and /i6200— 0.74, T p h — 4573 K, 
AEW qu i e t^l.728 A, and -F plagc /_P chrom = 4 have been used. 



Peg 



0.35 r 
0.30 : 
0.25 r 
0.20 : 
0.15 : 
0.5 
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0.7 



♦♦♦♦♦♦ 

0.8 0.9 



1 .0 



Fig. 7. Grids of solutions for II Peg. The filled circles rep- 
resent the values of spot temperature and area for the 
light-curve solutions, while the diamonds represent the so- 
lutions for the temperature curve. The hatched area, in 
each box, is the locus of the allowed solutions accounting 
for the data uncertainty. 



6. Discussion 

The solar paradigm can serve as a guide for the inter- 
pretation of the behaviors of stellar ARs at different at- 
mospheric layers. In the Sun, it is found that spots are 



always spatially associated with faculae, although the re- 
verse is not always true. Indeed, faculae generally form 
before a spot or a spot group appears and survive for some 
time after the disappearance of the spots. Moreover, high- 
latitude small faculae or plages not associated to spots are 
normally observed on the Sun. The ARs trace magnetic 
flux tube emersion, that is influenced by several effects, 
such as the Coriolis force, magnetic tension, and convec- 
tive motions. Therefore, studies of the topology, motions, 
and orientations of ARs at different levels (spots, plages) 
can provide information about the physical processes oc- 
curring below the solar/stellar surface. 

Some papers dealt with the behavior of facular 
to sunspot areas along the sola r activity cycl e (e.g., 
IChapman et al.lll997t lFoukal [T998). In particular. iFoukal 
(|1998l ) found that the area ratio of faculae to spots de- 
creases at increasing activity levels both using white-light 
faculae and Ca II K chromospheric plages and explains this 
as a consequence of the different dependencies of plage and 
spot lifetimes upon their emergent magnetic flux. Thus, 
the sub-photospheric field properties are believed to be 
more important to determine this ratio, rather than pho- 
tospheric field diffusion. Moreover, this shift toward dark 
photospheric structures for high activity levels can explain 
the high variation amplitudes observed in late-type stars 
more active than the Sun. 

Some authors found that for moderately active stars 
the variability at op tical wavelength is strongly influ- 
enced by faculae (e.g., Radick et al.lll998l : lMirtorabi et al 



2003). However, the facular contribution has been pro- 



posed also for very active stars by other authors to ac- 
count for the UV fl ux excess of a ctive stars compared 
to non-active ones ( Amado l2003h and to explain the 
blueing of the colors when the stars get fainter (e.g., 
Aarum Ulvas fc Engvoldlll999l : iMessina et al.ll2006l ). The 



presence of faculae/plages arou nd spots seems to be an 



ubiquitous phenomenon (e.g.. iFrasca et al.l I19971 I1998 



Catalano etail l2000l : iBiazzo et all I2006L l2007h and the 



non-detections could be due to contrast reasons. 

Despite the uncertainty in deriving the radii of the 
plages from our Ha variation curves, we would like to 
outline that the plage area is larger than that of the un- 
derlying spot for the smaller AR, while the reverse is true 
(the spot bigger than the plage) for the larger AR. This 
holds true bot h for A And a nd II Peg and is in line with 
the findings of lFoukall (|1998l ) for solar ARs. 

As regards the degree of spatial correlation of spots 
and plages in active stars, guidelines can be represented 



by dis k-integrated observation of the Sun. iCatalano et al 



(ll998l ) report on a strong rotational modulation of the so- 
lar irradiance in the C 11 A 1335 A chromospheric line from 
UARS SOLSTICE experiment. They ascribe to chromo- 
spheric plages this modulation and show that it is highly 
correlated with the sunspot number. The average position 
of spots appears to alternately lead and lag the centroid 
of plages by about 30-40°, at maximum, with a possible 
period of 270 days. Some indication of small longitude 
shifts between Ha plages and spots within this range has 
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been found for very active stars by Catalano et al.l ( 2000h . 
We find that the plage longitudes in II Peg are nearly the 
same as those of the underlying spots. For A And, instead, 
longitude shifts of about +17° and —16° for the two ARs, 
that we consider as marginally significant due to the data 
scatter and phase coverage, come out from the model so- 
lution. 

The investigation of possible dependencies of spot pa- 
rameters like temperature and filling factor on global stel- 
lar parameters like effective temperature, gravity, activity 
level (rotation rate, differential rotation, etc.) is of great 
importance to better understand the physical mechanisms 
at work on the formation and evolution of ARs. 

Recently, a correlation of temperature difference be- 
tween the quiet photosphere and spo ts, AT, with the e f- 
fective temperature has been found bv lBerdvuginal (|2005l ). 
She found that, on average, AT is larger for the hot- 
ter stars, with values of nearly 2000 K for stars with 
T GfT w 6000 K and falling down to about 200 K for M4 
stars. This behavior is displayed both by giant and main- 
sequence stars. 

It is interesting to investigate the role of the sur- 
face gravity on AT by selecting stars of nearly the same 
temperature. To further investigate this issue, we have 
also used the values of AT derived by us in Paper II 
for the three active giant/subgiant stars with the LDR 
method. In addition, we have also considered the values 
of AT— 1325 K and 1030 K obtained, w ith a typical error 
of ±150 K, by lO'Neal etal] (|200lLl2004h with the method 
of TiO bands for the two main sequence stars V833 Tau 
and EQ Vir, which have T e g similar to the star inves- 
tigated by us and \ogg ~ 4.5. The data of these seven 
active stars suggest an incr easing trend of AT with grav- 
ity, as already proposed bv lO'Neal et al. ( 19961 ). that can 
be explained by the balance of magnetic and gas pressure 
in the flux tubes of active regions. However, more data are 
needed to further investigate this correlation. 

7. Conclusion 

We have presented a study of the surface inhomogeneities 
at both photospheric and chromospheric levels based on 
a contemporaneous spectroscopic and photometric moni- 
toring of the two active RS CVn stars A Andromedae and 
II Pegasi. 

From the same data set of medium-resolution optical 
spectra, we have obtained information about the chromo- 
spheric and photospheric surface inhomogeneities by us- 
ing the Ha emission and the photospheric temperature 
(from line depth ratios), respectively. Additional infor- 
mation coming from the light curves, together with the 
temperature modulations, allowed us to disentangle the 
effects of starspot temperature and area and deduce these 
parameters in a unique way. A very important indication 
from this work is that the the starspots of A And are 
considerably smaller and warmer than those of II Peg, 
notwithstanding the nearly equal photospheric tempera- 
ture. At present we cannot say if this is due to the different 



gravity of the two active stars, being A And a giant with 
logy ~2.5 and II Peg a subgiant (log 3 —3.2), or if it is 
simply the effect of the higher activity level of II Peg com- 
pared to A And. More active stars with different gravity 
and activity level must be investigated to settle this point. 
However, by using the values of temperature difference be- 
tween photosphere and spot, AT, for other stars from our 
previous works and from the literature, we find an increas- 
ing trend of AT versus logg that could be explained by 
the magnetostatic equilibrium between gas and magnetic 
pressure. 

We have found, for all the stars observed, a tight 
anti-correlation between the Ha emission and the pho- 
tospheric temperature modulations, that indicates a close 
spatial association between photospheric spots and chro- 
mospheric plages. The largest longitude shifts between 
plages and spots of about ±16°, have been found for 
A And. Moreover, the area ratio of plages to spots seems 
to decrease when the spots get bigger. 

Furthermore, in II Peg a strong flare affecting both the 
Ha and Hei lines has been observed. The energy losses 
in these lines have been also evaluated. A possible flare 
event, with a much smaller energy budget, seems to have 
occurred in A And. In the present work we have shown 
the great power of a coordinated photometric and spec- 
troscopic monitoring of active stars for the study of the 
main properties of their active regions. 
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